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NASA MINI TRACK INTERFEROMETER 
R E F RA CT IO N CO RR ECTIONS 

P. E. Schmid 


ABSTRACT 

The NASA Minitrack interferometer system is an 
accurate, reliable and proven spacecraft tracking sys- 
tem. This report reviews the influence of the Earth's 
atmosphere upon Minitrack calibration and spacecraft 
tracking. It is believed that this review will aid the or- 
bit computation engineer in selecting those refraction 
correction procedures which will assure maximum orbit 
computation accuracy whenever an interferometer such 
as the Minitrack system is employed. 
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NASA MINITRACK INTKKKK HOM I TK H 
RKF HACTH )N COKRKC HONS 


P. I . Schmid 


SUMMARY 

This report examines the retraction effects of tlu* 
Karth's atmosphere on the NASA 130 Mil/ Minitrack 
interferometer during aircraft calibration and space- 
craft tracking. It is shown that onlv the troposphere 
(sen level to 30 km) affects aircraft calibration. It is 
also seen that the ionosphere (s5 km to 1000 km above 
Karth's surface) is the principal refraction bias source 
during spacecraft track. This results from the tropo- 
sphere eoi rec tion inherent in an interferometer angle 
calculation. 

The daytime ionosphere induced bias at 25 eleva- 
tion can be as large as 1 mr. This bias decreases to 
0.1 mr at M) elevation. Nighttime, when total ionosphere 
electron content decreases, generally results in onh 10 
of the daytime refraction. Possible exceptions have been 
cited for near equatorial stations where nighttime anomol- 
ous ionospheric behavior has caused Minitrack signal 
dropouts of telemetry and tracking data. This anomaly 
when present has been noted to generally peak at local 
midnight. It is clear that ionospheric refraction cor- 
rections are a necessary part of Minitrack data process- 
ing if daytime data below elevation angles of -43 is to be 
used for precision orbit computation. This is required 
to reduce the ionospheric refraction bias to a value com- 
parable to the Minitrack system resolution of 0.1 mr. 
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NASA Ml NIT RACK I NT I H l I : ROM FT FR 

hi i uachon corrictions 


l.u INI R( MH’Cl U )N 


This report reviews the influence ol the Farth’s atmosphere u|.m >n NASA 
Minitraek measurements during system calibration and spacecraft tracking. 
The need lor such a review stems from tlu* increased number ol precision or- 
bits now being generated based upon Minitraek data (lor example Reference 1), 
coupled with a lack ol formal documentation in the area of Minitraek refraction 
modeling. The basic considerations presented in this report apply equally well 
to any short baseline radio interferometer tracking sysleni. 


Section 2.0 of this report briefly compares the angle information provided 
by an interferometer such as Minitraek with that obtained from steerable an- 
tenna systems. NASA steerable or monopulse trucking systems include; 


1. The Cioddard Range and Range Rate System (Reference 2) 

2. The Unified S Band System (Reference 2) 


1). IMlse C- Band radar (Reference I) 


I. and the Applications Technology Satellite Range and Range Rate Sy stem 
( Relerenc e 5) 

The Karth’s atmosphere can for purposes of radiowave propagation analysis 
be considered to consist of three layers. These lay ers, as indicated by Figure 1, 
are a function of altitude above the earth's surface and include the: 


1. Troposphere 

0 

h 

30 km 

2. Free- space Region 

30 

h 

85 km 

3. Ionosphere 

85 

h •- 

1000 km 


Section 3.0 presents a brief description of the frequency independent tropospheric 
angle correction equations appropriate to steerable antenna spacecraft tracking. 
Section 1.0 then indicates why the interferometer, to a first order, inherently 
corrects for tropospheric refraction when the baseline length is measured in 
free space wavelengths at the spacecraft transmitter frequency. In Section 5 
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Figurr 1 - Re la ti onsh i p Between Troposphere and Ionosphere 


it is shown t h:it the elicit ol the one way Doppler shill upon Minitruek measure- 
ment does not significantly detract from this tropospheric angle eorreelion. 


Seri ion .">.() ol this report reviews the philosophy ol the present Minitraek 
aireral't calibration procedure indicating the significance ol radio and optical 
refraction. I ((nations are derived which properly model the atmosphere lor 
purposes ol calibration computations. Those particular equations are incorporated 
in the Minitraek star background photo plate reduction computer program used by 
the Physical Science Laboratory, New Mexico State I’niversity for the Goddard 
Space Flight (’enter. 


Finally, section (>.(), indicates what improvements can be anticipated by 
modeling the ionosphen . Ionospheric ray bending decreases in proportion to the 
square ol transmission Irequcney (Relerence t»). Thus the distortion introduced 
b\ the ionosphere witli an interferometer operating :tt S-Band frec|iieneies 
(nominal 2 (ill/) is reduced by a factor ol approximately 2 on relative to that at 
l.'hi MHz. By contrast, the tropospheric bending is frequency independent in the 
frequency range of 100 MHz to 20 GHz (Reference 7). 



I inure 2 indicates the* overall Minitrack data handling and aircraft calibra- 
tion interlace. With regard to Figure 2, tins report is primarily concerned with 


1. the optica) and radio retraction corrections applied t<> the ai re rail cali 
b ration data and, 

2. tile inherent troposphere corrections associated with data delivered to the 
orbit compulation user and ionospheiie refraction during spacer raft 

» racking. 


Currently aircraft calibration flights are performed approximately twice 
per year for each ol the NASA Mind rack stations. The site locations are pre- 
sented in f'in'u re .‘1. The Minitrack and the Goddard Range and Range Kate 
systems are the primary tracking systems associated with the NASA Satellite 
tracking and Data Acquisition Network (STADAN). 


It should be mentioned that the primary purpose ol the aircraft calibration 
is to determine the ur-lield phase characteristics of the Minitrack antenna 
array, The possibility ol obtaining the same information via residual ana ’.vs is 
from actual Minitrack satellite orbit calculations is currently under investiga- 
tion (Reference !)). \ regression analysis oi this type appears feasible il multiple 
station, short re passes are considered, fneertainties in tin earth gcopotential 
function tend to mask treking svstem biases when more than one earth revolu- 
tion is calculated (Reference 10), 


2.0 INTERFEROMETER VERSES STEERABLE ANTENNA AN' ’ LE 
MEASHREM ENT 


The Minitrack system is an interferometer tracking scheme (References 
11 and 12) which provides angular data in the form of direction cosines describ- 
ing the position of a spacecraft relative to a given tracking station. It is a 
passive system where the spacecraft transmits at a nominal 13(> Mil/. Hence 
no ranging information is available and at present no attempt is made to utilize 
the available one way Doppler frequency information for range rate calculations. 
The relationships between direction cosines and tor? conventional definitions of 
azimuth and elevation arc given in Figure -4. As will be shown, the Minitrack 
direction c ine calculation does not provide a measure of either signal eleva- 
tion angle of arrival or line-of-sight elevation to spacecraft but rather an ele- 
vation angle between the two which to a first order is the angle of arrival 
corrected for tropospheric refraction. Elevation angle of arrival is the eleva- 
tion angle of the incoming, essentially planar, wavefront. Line-of-sight eleva- 
tion is with reference to the straight line segment joining the interferometer 
array center and the spacecraft. 
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Figure 2-Overall Minitrack Data Handling 
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I hi> i in runt rust to thf « lev ation an^lv <d arrival measure provided l>\ 
tr il l. in- s\8tem> emplo\inn steel al»l« antennas te.i;., the NASA t allied S Hand, 


Station 

Latitude 

Longitude 

Geodetic 

Height 

(deg 

mm. 

sec ) 

(deg. 

mm. 

sec.) 

(meters! 

Fairbanks, Alaska 

64 

S2 

18 61 

212 

09 

40 t6 

189 

Goldstone Lake, California 

35 

19 

46 56 

243 

06 

00 86 

921 

East Grand f ork-*, Minnesota 

48 

01 

2118 

262 

69 

21 06 

249 

Fort Myers, Florida 

26 

32 

53 78 

278 

08 

04 60 

9 

Blossom Point, Maryland 

38 

25 

49 91 

282 

14 

49 3 / 

6 

St John’s, Newfoundland 

47 

44 

28 94 

307 

16 

46 n 

112 

Quito, Ecuador 

-00 

37 

20 55 

281 

26 

16 62 

3578 

Lima, Peru 

-11 

46 

34 86 

282 

60 

69 14 

34 

Antofagasta, Chile 

-23 

37 

14 07 

289 

43 

38 32 

516 

Santiago, Chile 

-33 

08 

56 23 

289 

19 

52 88 

681 

Winkfield, England 

51 

26 

45.43 

359 

18 

13.87 

07 

Johannesburg, South Africa 

-25 

53 

00 98 

27 

42 

28 49 

1365 

Tananarive, Madagascar 

-19 

00 

25 21 

47 

18 

00 46 

1361 


(Reference 8) Note Not all of the indicated 

bite' du- cu.rently active 

Geodetic Positions on MERCURY Datum (Oct 19b9). 


a 6 378 166m 1/f ■ 298.3 

•‘<|uatorial radius flattening factor 

where: 


a 

j - polar radius of earth ellipsoid 

Figure 3-Satellite Tracking and Data Acquisition Ne*'vork (STADAN Minitrack Stations 
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Figure* 4-Definition of Minitrack Coordinate* System 


Jz cos a = SIN AZ COS E 
nt: COS 0 -■ COS AZ COS E 


t ioddard Ran^'e and l!an"c 
lion alible ol arrival is ol 


Hale, and Pulse C- Hand radar systems). The eleva- 
course available direetlx from an azimuth-elevation 


antenna mount. 


It ean be* obtained I >\ wax ol straight forward coordinate transformations 
xvhen the an^le data is measured with X- h or polar antenna mounts (Reference 
13). 


3.0 TROmSPlIKKIC A N’t ILK CORRECTION APPL1 KD TO STEERABLE 
ANTENNA TRACKING SYSTEMS 

The usual eorreetion scheme applied to elevation angle data obtained from 
trackin'.', *;\ stems employing steerable antennas is: 


E 


R 0 “ L : E 


e„ - : e. 


E 

1 


(i 


( ) 


r. obse i*vcd elevation angle <>1 arrival 
I ;i correction ac counting lor atmospheric bending 


K t is associated with t ropospheric bending 
K i 1 ' assoc iated with ionospheric* I ending 


It should he pointed out that l lie* el loot of the* ionosphere’s presenc e is al- 
ways to increase the error angle E aljove that due to tlu* troposphere alone*. 
Since the bending within the ionosphere is in a direction opposite to that within 
llu* troposphere, one might guess that tlu* effects arc* pailialh compensating. 

I hat tiiis is not the* ease, will lie* shown in section G.O. However, an intuitive 
answer can lie* obtained with j’eferenee to Figure 5. By considering a fixed angle 
ol arrival, E , one can observe what the error, E, will be lor various spac e— 
c raft altitudes. 


For a given angle ol arrival, the ray path is determined and hence the* 
line-of-sight angle E will always be bc*tween the tracking station and a point on 
the ray trace corresponding to the spacecraft altitude, h . As is suggested by 
Figure 3, tlie error angle ( E E - E ) will never be less than the error due to 
the* troposphere alone. The error* actually has a maximum when the spacecraft 
is within the ionosphere and tends to that due to the troposphere alone as the 
spacec raft altitude reaches an altitude of several thousand kilometers above 
the* earth's surface. 


The correction represented by E ( (equation 1) is frequency independent in 
the range 100 to 20000 Mil/., The correction represented by 1E 1 (equation 1) 
varies inversely as the square of frequency. At frequencies greater than 1 Gil/, 
and at elevations greater than 3 the ionosphere accounts for no more than 2' ( of 
the* total refraction (Reference 6). At 130 MHz, however, the tropospheric and 
ionospheric refraction effects can be comparable depending upon spacecraft 
altitude. The tropospheric first order elevation correction for tracking space- 
craft at altitudes greater than 100 km and E () • 10 with a steerable antenna is 
given by: (References 14, 15, 10). 


AE, N s (10 f ) cot E 0 radians 


( 2 ) 


where: AL and E n are as defined previously and N , the surface refractivity , 

t U c S 

is a convenient measure of surface index of refraction. 
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E 0 » ANGLE OF ARRIVAL (FIXED FOR A GIVEN RAY PATH)* CONSTANT 
E * LINE-OF-SIGHT ELEVATION ( i.e. E, AND E 2 RESPECTIVELY) 

ERROR ANGLE * AE * E 0 - E 

V HH *Sl • 

NOTE THAT AE, >AE 2 AND AE -► CONSTANT AS 

Figure 5— Increase in Elevation Angle Error Due to Ionosphere 
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(3) 


n # = 1 ♦ N b (10 6 ) 


The surface refractivity, N , can be computed from: 


N 


77. 6 
T 



4810 e RH 

— r~ 


(•«) 


where: 

T tempt ature in degrees Kelvin 

P total atmospheric pressure (millibars) 

RH relative humidity 

e saturation vapor pressure (millibars) 

Equations (1) and (2) can also be written in terms of direction cosines. For 
frequencies above 1 GHz the ionospheric contribution to ray bending is usually 
considered negligible and equation (1) becomes: 


E = E 0 -N s (10~ 6 ) cot E 0 
or 

cos E = cos [E 0 - (10‘ 6 ) N s cot E Q ] = cos E 0 [1 + N s (10‘ 6 )] 


since N g (10 6 ) cot E Q « E Q for E 0 > 5° 


(5) 


( 0 ) 


but from Figure 4 

cos a = sin AZ cos E 
m = cos /3 = cos AZ cos E 
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and letting 


m„ 


represent angle of arrival direction cosines with , m repre- 


senting "best estimate" of line-oi- sight direction cosines, it is seen that: 


' 1 1 + N b (10- 6 )] i Q (7) 

and 

m = [1 +N s (10- fc )] m 0 (8) 


Summarizing then, it is seen that a steerable antenna (i.e. monopulse) track- 
ing system when corrected for first order tropospheric bending effects, results 
in the angle of arrival direction cosine adjustments indicated by equations (7) 
and (8). The next section will show that tins tropospheric correction is inherent 
in the Minitrack direction cosine calculation. 


1.0 TROIKXSPHERIC ANGLE CORRECTION INHERENT IN INTERFEROMETER 
M E AS l IRE M ENTS 


An ideal radio interferometer measurement will be defined as one where 
the incident signal is planar and operation is within a stabile (i.e. time in- 
variant over measurement period) spherically stratified exponential troposphere. 
These conditions are for all practical purposes met when radio energy is re- 
ceived from a spacecraft at altitudes in excess of 100 km above the Earth’s 
surface and when f requencies in excess of 1 GHz are employed thus obviating 
consideration of the ionosphere. At lower frequencies such as the 130 MHz 
transmission of Minitrack signals, the ionosphere will also affect the angle 
measurement. However, tliis ionospheric effect can be shown to produce a linear 
addition to elevation angle measurement which will appear as a residual error 
once the troposphere has been taken into account. The following derivation shows 
that the ideal interferometer provides a 'built in" first order correction for the 
tropospheric effect and that tins correction is the same as is normally applied 
to elevation angle data obtained from steerable antenna trackers (i.e. equation 2). 


d.l THE BASIC MI NIT RACK MEASUREMENT 

As suggested by Figure 6, the basic measuremer , of r orthogonally 
intersecting interferometer such as Minitrack, are the !U _tion cosines K, and 
m . More precisely, however, the measurement is a ph comparison, which 
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Figure 6 — S i rr» p I if i ed Diagram of Minitrack Interferometer Spacecraft Tracking 


lor each antenna pair, provides a measure of radio path length difference. This 
is indicated in figure (i where L Si and AS 2 in actual operation represent radio 
path length differences rather than geometric differences associated with oper- 
ation within free space. . The object of all elevation angle corrections to radio 
tracking data is to approximate the geometric line of sight elevation angle between 
tracking stations and spacecraft as closely as possible. This line of sight is 
indicated by S of Figure 0. 


The phase measurement associated with the two interferometer antenna 
pairs is indicated in Figure 7. 

With reference to Figure G, the direction cosines , m 0 of the angle of 
arrival are given by: 


! _ 
0 ~ 


c o s 


a o = 


M. \ 

1 s 


m 


o 


cos 





\ 


b 


(9) 

( 10 ) 
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WHERE a Q AND £ 0 ARE RESPECTIVE INTERFEROMETER 
ANGLES CORRESPONDING TO THE WAVEFRONT ANGL: 
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Figure 7- Interferometer Phase Measurement 


where and M, are respective whole and fractional cycles of phase shift 
measured by the Minitrack system, s , is the signal wavelength as measured 
in the refractive media in the immediate vicinity of the interferometer, and b 
is the baseline separation of each of the interferometer legs. 


4.2 MINITRACK ANGLE COMPUTATION 

The information contained in the raw tracking data is Mj and M 2 to some 
ambiguity. At Goddard ambiguities are resolved and the phase counts are con- 
verted to direction cosines using baseline separation in free space wavelengths, 
k , rather than the actual wavelength indicated in equations 9 and 10 (reference 
30). Thus the angles a.' , ,3' calculated using free space baseline separation are 
not a measure of angle of arrival but some "pseudo angle" which is linked to the 
angle of arrival by the surface index of refraction. That is: 


k 


s 



n 

s 


( 11 ) 


12 


and hence 


cos I ' 


n, cos . 0 - [1 ♦ N, (IQ 6 )] i Q 


( 12 ) 


similarly 


m* - [1 ♦ N s ( 10 6 >] m 0 


( 13 ) 


where ' and m' are calculated direction cosines, and , and m 0 are angle of arrival 
direction cosines. However, equations 12 and 13 are seen to be exactly the re- 
sults presented previously by equations 7 and 8 where a first order tropospheric 
correction had been applied to a monopulse (steerable antenna) tracking system. 

Thus the Minitrack angle calculation has a "built in" tropospheric refraction 
correction which for E 0 > 10 accounts for at least 98'7 of the tropospheric 
bending. If no atmospheric corrections are applied to Minitrack data it can be 
expected that the primary refraction effect of the Earth's atmosphere will be 
linked to the ionosphere. 


5.0 OPTICAL AND RADIO RK FRACTION DURING MINITRACK CALIBRATION 

At present the Minitrack system is calibrated by means of aircraft which 
carry both optical and radio beacons. The optical beacon is photographed against 
a star background which permits a precise means for determining the line of 
sight between station and aircraft (References 17 and 18). 

The primary purpose of the aircraft flyby is to determine the far field 
phase characteristics of the Minitrack antenna array. Since Minitrack antenna 
spacing corresponding to each of the intersecting baselines is on the order of 
50 wavelengths, the far field should be probed at slant ranges of approximately 
10 km or greater. The aircraft also provides a "zero set" phase measurement 
which accounts for inequalities in cable electrical length between respective 
antenna arrays and the Minitrack receiver. These cables are filled with 
moisture free Nitrogen and maintained at an approximate gauge pressure of 5 
pounds per square inch. In addition to the aircraft flyby, internal Minitrack 
system calibrations are performed to account for electronic equipment phase 
delays. 
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Figure 8-0ptical Refraction Geometry 


5.1 CALIBRATION AND OPTICAL UK FRACTION 

As indicated in Figure 8, the calibration aircraft is photographed against a 
star background. The line- of- sight from the center of the Minitrack orthogonal 
pair of interferometers to a given star is then accurately determined since the 
star field has been well established by astronomers. However, the line of sight 
to the aircraft is different than the line-of- sight to the star by a small angle < . 
This angle, t , must be added to the line- of- sight angle to the star to provide an 
optical line-ol-sight to the aircraft which is then corrected for radio refraction 
to permit a calculation of what the interferometer should read at 136 MHz. An 
expression for the angle, < , has been rigorously derived in Reference 19. 

As will be shown in section G.O it is also given to a good approximation by: 

e = ( n 0 - 2 + — \ cot E (radians) (14) 

\ n o / 


14 


wlu*re: 


E elevation angle to the star 

n average optical index oi refraction Ijctwecn Minitrack station and 
aircraft altitudes 


1 



f H A 

n(h) dh 

1 


n (l surface optical index of refraction 

ll A altitude oi aircraft alxjve sea level 

H altitude of tracking station above sea level 

which for an exponential spherically stratified troposphere 
(i.e. n(h) 1 + N 0 (l(T fi ) e" vh ) leads to: 


N 0 (10" 6 ) cot E [1 - e' V “ A " #) ] 

= - » 7 > 


(15) 


The validity of the exponential model for purposes of radio refraction cor- 
rection has been well established (Reference 20, 14, 7). 

With regard to equation (15) the optical surface refraetivity N n can be 
calculated from: 


N 


77.6 p 
T 


(16) 


where 

p total atmospheric pressure (mb) 
T air temperature (°K) and 
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The deca\ constant, , for the optical Index of ref ruction within the I ailh's 
atmosphere is approximated 0.1 km' 1 (deference 19). 

I he next step is to determine what the Minitrack direction cosine measure- 
ments should he, based upon the line-of-sight angle* to the aircraft (K ol 
Figure 8). Lot cos , cos be the direction cosines to the star coincident with 
the flashing light and cos cos ' be the desired line- of-* sight direction 
cosines to the aircraft. Then: 


cos i ' cos i 


cos a 

COS I 


and 


(17) 


I /cos ' 

cos = cos 


'OS 


( 18 ) 


from previous definitions 


cos i sin AZ cos E 


cos - cos AZ cos E 


cos i r s i n 


AZ cos (E + «:) = sin AZ (cos E - « sin E) 


cos ' - cos AZ cos (E + t) - cos AZ (cos E - t sin E) 


(19) 

( 20 ) 


Thus the ratios cos / /cos i and cos '/cos (equations 17 and 18 respectively) 
are given by: 


cos a 1 cos 
cos I cos 


1 - € tan E = (1 - K) 


(M) 
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IH combining < 1 .“») and rJl), the correction factor (1 - k) becomes: 


and 


1 - K 


1 - 


-vh.-h.). 

N n (10 # ) [i .(■ A ■ 




(22) 


cos i ' = cos 1 (1 - K) 


(23) 


cos ' cos ( 1 - K) 


( 21 ) 


cos and cos being direction cosines to a star coincident with the flashing light 
and cos , cos ' direction cosines to the aircraft. 

Kquation 22 is what is currenth employed in the calibration calculation at 
the Physical Science Laboratory, New Mexico State I nivcrsitv. 

The next section will show that the expected radio phase meter readings M 
and ‘ . . are given by: 


and 



b cos i (1 - K ) n 



r ad i an s 


(23) 



lb cos 


(1 - K) n] 



r ad i a ns 


( 20 ) 


where: 


b baseline length, cos i,cos are direction cosines corresponding to 
line- of- sight to star occuiated by flashing light 

(1 - K) correction factor to obtain line- of- sight to aircraft 
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/» M, 


n - 


Ik - H I 

A * •'h 


n(h) ilh * average* radio index ot ref met ion brtwi’cn 
t r ark inn stat ion and n i r c r a f t a 1 t i t titles 


live- snare wavelength nl received frequency 


7, .2 C’ALI BUATION AND UADIO KF FRACTION 

During the aireralt calibration flyb> at :t given site, radio phase data is 
recorded and he tlilierence between the observed phase data and that "expected" 
is sent to the Goddard Space Flight Center as a correction to be applied to sub- 
sequent raw Minitrack phase data. During this calibration the "expected" Mini- 
track radio phase measurement is based upon the precise line oi sight to air- 
craft adjusted for prevailing radio refraction conditions. The previous section 
discussed the conversion ot stai line-of- sight to aircraft line-of-sight. This 
section will now derive the radio refraction correction fa 'tor. 

During calibration the ionosphere plays no pail since the aircraft is neces- 
sarily well below the ionosphere's 85 km lower altitude limit. A most important 
point is that the radio phase readings measured during calibration arc compared 
with precise line of sight data altered slightly to take into account the tropospheric 
r dio refraction. The calibration is not a means for later adjusting raw space- 
craft tracking data to either angle of arrival or line of sight measurements. The 
calibration accounts for on site biases associated with the Minitrack antenna 
system and electronics and does not provide a refraction correction for subse- 
quent raw data recordings. 

In order to calculate what the interferometer should read for a given line 
of sight geometry it is necessary to calculate the electrical path difference 
corresponding to the geometric path difference Sj - S, (Figure 5)). It will be 
shown that, for the aircraft geometry, the corresponding electrical path length 
is given by nJSj - S 2 ), where n is the average radio index of refraction between 
tracking site and aircraft altitudes. 

With reference to Figure 9, in the absence of an atmosphere, the free space 
interferometer phase measurement would be given by: 


2tt(S. - S ) 

A <? Q - (radians) 


( 27 ) 
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H. = APPROXIMATE AIRCRAFT HEIGHT (NEGLECTING 
EARTH CURVATURE) 


h j 1 station altitude 


NAVA Oil •- 

V *. K »«A « * >*• • ASA • 

ANAS. .. OAM , v, 

• * A I**"* HQ* S'’ 


Figure 9-Aircrcft Calibration Geometry 


where 


Sj geometric distance from phase center of antenna at one end of base- 
line to aircraft 

geometric distance from phase center of antenna at other end of 
baseline to aircraft 

free space wavelength c/f 0 and c is the speed of light in vacuum, 
the frequency of the received radio energy. 


The frequency, f 0 , will differ from the transmitted frequency b\ the one- 
way Doppler shift f 0 (r/c) where r is the radial velocity of the transmit beacon 
relative to the station. At 13G MHz this is entirely negligible relative to the 
expected uncertainty in measured surface index of refraction of ±1 N unit, that 
is, the ratio of AN (*1) to N s (~ 300) is always large compared to f /c . The 
physics linking surface refractivity to wet and dry bulb temperature measure- 
ments and statements of expected accuracy are covered in detail in Reference 7. 
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Again with reference to figure 9 the aetr.al interferometer phase readings 
are given by: 


<* ; - — — ^(s) dSj - J*n(s) dS^ radians 


As pointed out by Bean and Thayer of NBS (Reference 16), equation (2M) can 
be integrated assuming flat-earth conditions for elevation angles E 0 > 10°, slant 
range r < < n 0 sin E Q (a 0 earth radius), and for baselines of a few hundred 
meters or less. 

These conditions are applicable to the aircraft flyby and hence, again with 
reference to Figure 9, equation (28) can be written as: 


a- f H A ~ r H A 

r P = I n(h) dh I n(h) dh (radians) ( 29 ) 

\> sln >t J„ \ sin >2 J„ 


and by the law of sines 


sin 7 . s i n y sin 7 , sin y 


equation (29) becomes 


_ 2 ”( S , - S 2>/ 1 


H, - H 

v A s / * 


A 2tt(S - S ) 


pj - (radians) 
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With reference to figure 9, Sj - s,, at aircraft slant ranges, is given by: 


S j - S , = b cos . 


for one Minitrack baseline d , for the other orthogonal Minitrack baseline 
' and relating these to the star direction cosines, i and , equation (30) 
can lie written as 


2rr 

P« 8 ’’ cos a (1 - K) n — — (radians) 

o 

(31) 

Aqj . = lb cos (1 - K) n 1 -I— (ra d inns) 

"o 


where the factor ( 1 — k) was defined previously by equation (22). 

In equation (31) Acp a and Atp* represent the expected Minitrack phase read- 
ings during the aircraft flyby when the flashing light appears coincident with a 
star whose line-of-sight direction cosines are cos a and cos ■ respectively. 

The average radio index n is calculated using the NBS exponential profile 
n(h) = 1 + N s (10 6 ) e kh . This profile when integrated and averaged results 
in: 


n 


1 + N s (1C)" 6 ) (1 - e k ^ A 

k ^A- H s) 


(32) 
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the* dcca\ factor, k, is calculated with an empirical relation formulated In NHS 
based upon a lit to 888 sets of data (Reference 14). 


\ 



s 


where N represents the change in 
above the tracking station. 


relractivity in the vertical direction 1 km 


A N — 7 3 ' i , 0 . 0 0 5 5 7 7 n s 


(Refer one e 14) 


(lid) 


Finally the surface relractivity is calculated from surface meteorological 
measurements using the relation: 


N 


/ 77.6 \ / p 481 Of (RH) 
\ T I \ T 


(34) 


where n = 1 + N s (10 -6 ) surface radio index of refraction and 

T air temperature ( C K) 

P total air pressure (millibars) 

€ saturation vapor pressure (millibars) 

RI1 relative humidity (Reference 21) 

The optical relractivity, N 0 , as previously indicated by equation (10) is 
equation (.34) without the water vapor term. 

Equations (31) through (34) are used by New Mexico State University in their* 
reduction of Minitrack aircraft calibration data for the Goddard Space Flight 
Center. 
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5.3 ANC'» LF OF AKK1VAL CORKFCTION KQUATIONS 

As a matter of interest, the foregoing principles can now be applied to de- 
rive the flat earth approximation for the angular difference between elevation 
angle of arrival and line- of- sight for both optical and radio refraction. The only 
difference between the optical and radio cases is that n for the radio refraction 
includes the effect ol water vapor within the earth’s atmosphere. 

For an interferometer such as Minitrack which has intersecting orthogonal 
baselines running Fast-West and North- South the phase readings are, as shown 
previously, given by: 


and 


2rrb cos E — . 

n sm A 

* \> 



27rb cos E . 

n cos A 


(55) 


Here, E, represents the line- of- sight elevation angle to the source and, A, 
the azimuth angle. But the phase measurement is also represented by the phase 
delay in terms of the number of wavelengths of delay within the media surround- 
ing the antenna. The effective wavelength within the media in the vicinity of the 
antenna is \ 0 /n s . In terms of the elevation angle of arrival, E n , the following 
equations can thus be written: 


A<p a 


2 d) cos E 

\> n s 


0 


sin A 



2vb cos E 0 

X Q fl S 


cos A 


( 30 ) 


Combining (35) and (30) 


n cos E = n cos E 0 


(37) 
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This equation links elevation angle of arrival, K , to line-of- sight angle, K. Dut 
E n E E, where E is a small perturbation and (37) ean be expanded to ob- 
tain t lie approximation: 


E - 



cot E 


(38) 


Here E is the elevation angle- of- arrival error clue to the atmosphere. This 
result was also derived by R. H. Paul (equation 19 oi Reference 1 G ) using a 
somewhat different approach. Equation (38) can also be written in terms of 
refractivity, N, since* 


and 


n = 1 + N (10“ 6 ) 


n s = 1 + N s (1C)" 6 ) 


In all cases the deviation from unity is It)" 3 or less thus equation (38) can 
be written as 


AE = (N s - N) cot E (10" 6 ) 


(39) 


Where as before, N , represents the radio surface refractivity varying from 
200 to 400 depending upon geographical location and local meteorological condi- 
tions and N is the average refractivity between the station altitude and the alti- 
tude of the signal source. If the source is a spacecraft than N includes the effect 
of the ionosphere where the absolute sign of the refractivity term is negative 
(Reference 0). Equation (39) thus also indicates how the presence of the ionosphere 
increases the total ray bending as suggested previously by Figure 5. / the 
altitude of the spacecraft increases much above 5000 km, or for all spacecraft 
operating at 2 GHz or higher such that the ionosphere is no longer a significant 
influence, N tends to zero and equation (39) reduces to: 
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) iot E 


AE = N ( 1 0" 6 

1 ' 


and to a first order: 


cot E - cot E 0 


or 


AE = N s (10“ 6 ) cot E 0 (radians) (40) 


This result, presented earlier as equation(2), isthe usual tropospheric bendingcor- 
rection applied to ele\ ttion angle spacecraft tracking data for elevation angles 
above 10°. The same result is used for optical correction to star-light angle of 
arrival except N 0 , the optical surface refractivity , is used instead ol N s , the 
radio- surface refractivity. 

Finally, equations (38) and (40) can be used to derive equation (14) presented 
earlier for the angular difference, <- , between the line-of-sight to a star coinci- 
dent with a flashing light on the Minitrack calibration aircraft and the line-of- 
sight to the aircraft. 

With reference to Figure 8, the angle, e , is the difference between the total 
star-light bending, N o (10‘ 6 ) cot E, and the bending of the flashing light ray path 
relative to the line-of-sight to the aircraft, (1 -n/n 0 ) cot E or 


or 


€ - 



/ 




cot E 



(radians) 


(41) 


Equation (41) was presented earlier as equation (14) and presents a good 
approximation of elevation angle distortion for elevation angle, E 0 , greater than 
10° and slant ranges to the spacecraft much less than a 0 sin E 0 , where a Q is 
the earth radius. 
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0.0 III 1(1 ( )1 IONOSPII I HI ON MINITRACK MEASUREMENT 


The elevation angle bin: present in Minitraek measurements below elevation 
angles ol 10 can be largely attributed to residual error alter nominal ionospheric 
eorreelions have been applied. Unlike the troposphere, the ionosphere is much 
less predic table and as menti ic'd previous ly , ionospheric refraction effects var\ 
inversely as the square of the radio frequency. One of the major sources of 
ionospherie modeling error has to do with the 1 fact that the radio signal often 
traverse's portions of the ionosphere several thousand kilometers from the track- 
ing site 1 , a situation whole the usual spherically stratified Chapman model does 
not strictly apply . One means ol minimizing the ionospheric refraction error 
is to use only data where the signal is traversing local night regions. Recent 
studies, however, indic ate ionospheric propagation anomalies peaking at local 
midnight at t he* South American Minitrack site's (Santiago, Lima, and Quito) 
(Reference 22). Those equatorial region anomalies appear to predominate in 
the 1 Hast- West baseline whic h, as stated in Reference 22 is consistent with 
existing theories of ion plasma density alignment North to South along the 
Earth’s magnetic field lines. Study of the phenomena and the effect on NASA 
120 MHz telemetry , tracking and command signals is continuing (References 
22 and 24). 


<1.1 THE INTERFEROMETER REFRACTION ERROR EQUATION 


The inherent tropospheric correction to interferometer data was discussed 
in section 1.0. It was shown that the elevation angle calculated is given by: 


E' = E. - N (10“«) cot E,, 

Os ' 0 


(42) 


E' r E 0 - N s (10' 6 ) cot E 


for E Q > 10 where: 

E 0 elevation angle of arrival 

E' interferometer calculated angle 

E line-of-sight elevation angle 
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It was also shown in section 5.3 (equation d*) that the line-of-sight elevation 
angle, I , is linked to angle of arrival elevation, K (1 , by: 


E = E 


o 



Mil) 


Thus the total interferometer error is found by combining M2) and ( Id) to 
obtain: 


In t e r f e r ome ter Error = E' - E ~ - N ( 1 0 6 ) co t E 


Ml) 


which corresponds to Bean and Thayer's equation 2(i in Reference 10. liquation 
( I i) is a good approximation for E > 10°, interferometer baselines of a few 
hundred meters or less (such as Minitrack), and slant range to spacecraft, 
r < < .i Q sin E, where a 0 is the earth radius. Equation (44) is useful in 
qualitatively describing the effect of the ionosphere on Minitrack data when the 
spacecraft is within the ionosphere (100 < altitude < 1000 km). At ldG MHz the 
contribution to N’ due to the troposphere is negligible compared to the ionospheric 
contribution. This can be shown as follows. According to equation (d2) for alti- 
tudes significantly above the troposphere scale height (~7 km) the tropospheric 
average refractivity is given by 


N (10 -6 ) 

N = — (4( 

1 kH 

where: 

H spacecraft height above the tracking station 

k decay factor in exponential model, 1/k 7 km. 

For H 200 km, N 12 which is small when compared to the ionospheric 
cont ribution. 


The ionospheric refractivity is given by: 


-40.3 N (h) 

(1 0* ) ( Rc f o r cnci 1 6 ) (46) 

f a 

where: 

N electron den. trons/meter 3 ) 

f 9 

f frequency (Hz) 

For a daytime peak density at the height of maximum ionization, N takes 
on values of typically 0,5 ( 1 0 1 2 ) electrons/meter 3 and at 136 MHz this cor- 
responds to a maximum ionospheric refractivitj on the order of -2000, This 
value will taper off by a factor of 10 at the lower edge of the ionosphere (II 85 
km). However, the average contribution will be significantly greater than the 
tropospheric contribution N ( 12. An interesting observation is that, again by 
equation 1 1, and noting that the average refractivity due to the ionosphere is 
negative, the ray bending will be in the same direction as that due to the tropo- 
sphere. With the interferometer, of course, the tropospheric bending is aiieady 
inherently corrected for by the usual angle calculation discussed previously in 
section 4.0. 

The "flat-earth" model (equation 44) is useful for qualitative ionospheric 
refraction discussion only since the range of validity includes the constraint 
r < < a sin E. However a more satisfactory equation which includes Earth 
curvature has been presented by Be. n and Thayer of the National Bureau of 
Standards where the interferometer error is given by: 



r > > b b < 50 km 
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where: 

K line-of-sight elevation angle 

N(h) vertical rofractivit> profile Including both troposphere and ionosphe re 
i slant range station to spacecraft 

II altitude of spacec raft above station 

a 0 earth radius 

b interferometer baseline (b 120 meters lor Minitrack). 

For the "flat earth" approximation r sin E II and i ,<r < a 0 sin K, equation 
(17) reduces to -N(10“ f ’ ) cot E, presented previously as equation (11). 

Equation 17 can be derived from Snell's law for a spherically stratified 
media. The result can also be obtained by ray tracing, through the atmosphere 
(again using Snell's law) to determine the difference between angle of arrival 
and line-of-sight. This can then be linked to the interferometer error by sub- 
tracting out the first order tropospheric correction of N h (10” f> ) cot E. 


<>.2 MiNITRAUK IONOSPHERIC RE FRACTION 

This section presents the Minitrack bias error which can be expected if 
ionospheric effects are neglected. A nominal correction based upon National 
Bureau of Standards (NBS) Monthly worldwide peak electron density predictions 
and some mathematical distribution model (such as a Chapman profile) often 
reduces the error by only 50', . As will be discussed shortly the greatest un- 
certainty is related to the F, layer vertical sounding (i.c. "zero") Maximum 
Usable Frequency (MUF). The peak electron density and hence ref inactivity is 
proportional to the square of the F, (zero) MUF. 

As indicated by Figure 5 and equation 47. the Minitrack refraction error is 
separable into tropospheric and ionospheric components. At frequencies such as 150 
MHz , where ionospheric bending is appreciable the tropospheric contribution in 
equation 47 is generally negligible if elevation anglesabove 10 and spacecraft altitudes 
above 100 km are considered. During aircraft calibration the reverse is true and only 
the troposphere is a consideration. Bv ray tracing (Reference 25) through the 
ionosphere and troposphere separately and then comparing with a trace through 
ionosphere plus troposphere, the validity of superposition of bonding has lieen 
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verified. I hr difference between total bending due t<> the troposphere alone and 
N (id”' ) rot I comprises the tro|x>sphcric contribution to the error in equa- 
tion IT. The ionospheric caused elevation on or in equation 17 ran be obtained 
di recti' lrom a ray trace through the ionosphere alone. 

I able I is an example <>| the Minitruck refraction error magnitude associated 
with a track ot the 2000 km altitude geodetic satellite (GKOS- A) which was 
launched on 0 Noven. her 1905. It should be emphasized that exact refraction 
values will van from site to site and depend upon the time of day as well 
as time of year. For the table a peak ionospheric density of 0.8 ( 1 0 1 2 ) electrons/ 
meter 3 was assumed for daytime and 1 0 1 1 electrons/meter 3 for nighttime. 

TABLE 1 

Minitrack Refraction Error 
(No Corrections Applied) 


Elevation 

Angle 

Tropospheric 
Bias (mr) 

Ionosphe ric 
Bias (mr) 

(Degrees) 


n 

Day 

Night 

10 

-0.0G 

2.25 

0.30 

15 

-0.02 

1 .05 

0.20 

20 

-0.01 

1.25 

0.15 

30 

-0.005 

0.80 

0.10 

40 

-0.002 

0.50 

0.05 

GO 

-0.001 

0.25 

0.00 

80 

— -0.0005 

0.10 

0.00 


Satellite Altitude 2000 km 


i 


In Table I the term "bias" refers to that value which must be subtracted 
from the measured elevation to obtain "true" elevation. The minus sign on the 
troposphere bias simply indicates the extent to which the inherent N s (10~ 6 ) 
cot K correction of the ; liter ferometer overcompensates for the tropospheric 
bending. Also the ionospheric refraction is a maximum for satellites at an 
altitude of approximately 500 km. For an altitude of 500 km the ionospheric 
error values of Table 1 will be increased by approximately 25 r /c at 10 r , 50', at 
20 and 70' at 30". The ionoconeric bias tabulated in Table I is also plotted in 
Figure 10. The parameters used to obtain the results of Table I and Figure 10 
via a ray trace are as follows: 


. 



ELEVATION BIAS (MILLIRADIANS) 


r 


MiNiTRACK IONOSPHERIC PARAMETERS 

(JANUARY 1966 PREDICTION) 



Figure 10-Minitrack Elevation Angle of Arrival Bias Due to 
Ionosphere if no Correction Applied 


TROPOSPHERE: 


NBS exponential model 


N(h) = N s e" kh 

N = 3S0 k = 0.16 km" 1 


n (b) = 1 + N(h) (10“ 6 ) 


0 < h < 50 km 


(18) 
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ionc >s i»n i:m 


Chapman Model 

-40.3 N (h) 

N(h) = 10 6 

f 2 


f 136 MH/ 


n(h) = 1 + N(h) lCf 6 h > 50 km 


N 


(h) = N e* 1/2(1 ’ z ' e " Z) Reft* 


rcncc 
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(49) 


N maximum electron density 0.8 (10 12 ) electron meter 3 representing 
nominal '•'vtime value, 



h m 300 km representing a typical daytime height of maximum 
ionization 

H = scale height = 83 km 
h - 2000 km (GEOS- A altitude) 


For a more precise determination of refraction effects all of the foregoing 
parameters must he estimated in view of particular station geographical loca- 
tions, time of year, time of day, phase oi sun spot cycle, prevailing meteorologi- 
cal conditions, and so on. It should be pointed out that at present, except for 
special tracking tests, only Minitrack data above elevation angles of 80° is used 
for orbit computation. 

Figures 11 and 12 are plots of elevation angle residuals (observed elevation 
minus orbit calculated elevation) when no refraction corrections were applied 
to a 7 day stretch of GEOS- A Minitrack data (31 December 1965 to 6 January 
1966). It is felt that the trended portion of these curves can be primarily at- 
tributed to refraction within the ionosphere. 
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ELEVATION RESIDUAL 
(MILLIRADIANS) 


I 5 



NEWFOUNDLAND 
I5H < t< I9H 

LOCAL time 


FORKS (NORTH DAKOTA) 
I5H < t < I9H 
LOCAL TIME 


MOJAVE (CALIFORNIA) 
I5H < t< I9H 
LOCAL TIME 


Figure 11-Minitrack E levation Angle Orbital Residuals 
GEOS A (No Ionospheric Correction) 
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ELEVATION RESIDUAL 
(MILLIRADIANS) 


FORT MYERS (FLORIDA) 
I5H < t < I9H 
LOCAL TIME 




BLOSSOM POINT (MARYLAND) 
I5H < t< I9H 
LOCAL TIME 



COLLEGE (ALASKA) 
I7H < t< 20H 
LOCAL TIME 


NASA i,S» T%u*. 

M N THA l ' H> ANA.IM. N 

■ •»AN( M 551 OAT* i 

B* *■ ~ «i ►LOT NO ’ ■ 4/ 


Figure 1 2--Miriitrack Elevation Angle Orbital Residuals 
Gt OS A (No Ionospheric Correction) 


Finalh , a most important point is that the* tropospheric error indicate*.! in 
Table I would be the only significant atmospheric refraction bias for an inter- 
ferometer operating at frequencies above 2 GHz. For such a situation and for 
satellites where the slant range approaches the magnitude of the Faith radius 
multiplied by sin E, the interferometer refraction bias is given by the term 
associated with l/n 0 sin* F of equation 17. This is the expected interferometer 
refraction bias presented by the National Academy of Sciences (Reference 27) 
for interferometers with a baseline less than 3 km and operating at high enough 
frequencies such that the ionosphere is of no concern. The results in Reference 
27 check those for the tropospheric refraction presented in Table 1, the latter 
having been obtained from ray trace considerations. 

One might wonder why 130 MHz is used instead of say 2GHz for the Mini- 
track frequency to eliminate ionospheric refraction. The primary reason is 
historical since in l ( jr>6, when Minitrack was first implemented, efficient S-Band 
spacecraft qualified hardware was not available. However the use of 13(5 MHz 
still has certain advantages over S-Band such as ease of signal acquisition and 
ambiguity resolution. That is, the transfer of radio energy from a spacecraft 
using a low gain antenna to a given ground effective antenna aperture is indepen- 
dent of operating frequency (Reference 28). However the respective beam and 
lobe widths vary inversely as the operating frequency. The antenna patterns 
corresponding to the Minitrack fine resolution array as presented in Reference 
29 are given in Figure 13. Note that there are two such arrays for each of two 
orthogonal baselines. The complete Minitrack antenna layout as given in Refer- 
ence 30 is shown in Figure 14. 

0.3 IONOSPHERE MODEL PARAMETERS 

The Chapman model presented in section 0.2 as equation 49 requires three 
input parameters - namely, the maximum electron density, N , the height of the 
maximum, h , and the scale height, H. Of these three the most critical param- 
eters is the maximum density. The scale height has been formulated in terms 
of h (Reference 31). 


H = 1.66 [30 + 0.2 (h - 200)] 


where H and h are in km. 


The height of the maximum, if no other information is available, can be 
taken as h 300 km for daytime and h 250 km for nighttime (Reference 32). 
The overall refraction is not a critical function of h and H. 









Figure 14-Relative Location and Function of Antennas in a 
Mimtrack Antenna Field. 


The following indicates a means for estimating the peak electron density of 
the ionosphere, utilizing the Institute for Telecommunication Sciences-Environ- 
mental Science Services Administration (formerly Central Radio Propagation 
Laboratory-N3Sj monthly publication, "Ionospheric Predictions," in conjunction 
with a world-wide "gyro frequency" chart. 

The publication, "Ionospheric Predictions" presents F 2 (zero) MUF (maxi- 
mum usable frequency) world-wide charts (reference 33) for 3 months in advance. 
The charts present monthly mean values in 2-hour increments since the ioniza- 
tion variations are primarily diurnal during any given month. The "zero" in F 2 
(zero) refers to zero ground "skip distance," i.e. vertical incidence ionosphere 
sounding. It can be shown (see for example, Reference 34) that for vertical inci- 
dence the critical freouency (highest frequency signal which will be returned to 
earth) is a function only of peak electron density and magnitude of the earth's 
magnetic field in the region of reflection. The direction of the earth's magnetic 
field does not influence the vertical "critical frequency." Because of the 
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presence of the earth's magnetic field, 3 discrete critical frequencies exist - 
namely, the ordinary critical frequency (this frequency would be the only 
critical frequency fn the absence of the earth's magnetic field) and 2 so-called 
"extra-ordinary" critical frequencies. The F, (zero) Ml F charts are for the 
highest critical frequency or the upper extra-ordinary frequency. This upper 
frequency results in the more predominant reflection since signals at the lower 
frequencies encounter appreciable attenuation. 


The peak electron density is obtained from ihc ordinary critical frequency 


as: 


f 2 

N = — 
e 80 


(51) 


where: 


f ordinary critical frequency (Hz) 

N peak electron density (electron/meter 3 ) 

e 

The F (zero) MUF chart is for the upper extra-ordinary frequency (e.g. 
Figure 15) and is related to the ordinary critical frequency as follows: 


where 



(52) 


f upper extra-ordinary critical frequency as obtained from F 2 
1 (zero) MUF chart. 

f c ordinary critical frequency required to calculate maximum N p 
(Equation 51). 

f "gyro frequency" which couples the effect of the earth's magnetic 
field into the ionosphere index of refraction. (Figure 15) 

1 "7 r 

= — — B MHz (Reference 34) 

2v 

B = earth's magnetic field (Gau. s) 
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Since 17/1 I' in most cases then: 

h c 


f 

x 


1 



or the ordinary critical frequency is simply: 


f 


r 




(53) 


The predicted and observed F 2 (zero) MUF values seldom differ by more 
than 1 or 2 MHz with the larger differences occur ling at the beginning or ending 
of a given month. Unfortunately the maximum electron density vaiies as the 
square of the critical frequency (Equation 51) so for a value of f c 5 MHz under 
worst case conditions (error of 2 MHz) the electron density could be off by a fac- 
tor of 2. If the predicted critical frequency is off by only 1 MC then the maxi- 
mum electron density in this case would be correct to within 30', . During mid- 
month, the F 2 (zero) MUF predicted values are accurate to within 80 to 9 0 c /< . 

The maximum electron density is given by: 




electrons me t e r 3 


(54) 


f frequency for appropriate geographical location, month, and time of 
1 day as obtained from an F 2 (zero) MUF chart such as Figure 15. 

f gyro frequency appropriate to a particular tracking site (Figure 10). 
Note that the gyro frequency is not a function of time. 

N m estimated maximum electron density in the region above a particular 
tracking station (eiectrons/meter 3 ). 
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JANUARY ItfW) UT 07 I GMT I 

longitude 

90° 120° 150° 180° «0° 120° 90° 60° 

EAST WEST 


WEST EAST 


90° 120° 150° 180° 150° 120° 90° 60° 

EAST WEST 

PREDICTED MEDIAN MUF (ZERO) F2 (MHi) 
(PLOT OF f. ) 


WEST EAST 


(REFERENCE 33) 


ero) MUF Prediction (Reference 33) 


The maximum electron density, N.. calculated by equation 53 is used in the 
Chapman model (equation 49) as a scale factor for the vertical electron density 
profile. To obtain best results the profile should lx? calculated for the region 
where the radio energy intercepts the region of maximum density. Generally, 
however, only the region above the tracking site is considered and spherical 
symmetry is assumed. The error thus introduced depends on satellite altitude 
and the relative steepness of contours such as given by Figure 15. 


7.0 CONCLUSIONS 


The atmosphere, for purposes of radiowave propagation, can be separated 
into two parts - namely, the troposphere (sea level to 30 km) and the ionosphere 
(85 km to 1000 km). Troposphere ref: iction effects are frequency independent 
in the range 100 MHz to 20,000 MHz. Ionospheric refraction is frequency 
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SOUTH north 


EAST WEST WEST EAST 

ISO” 180° 150° 120° 90° 60° 30° 0 ° 30° 60° 



igure 1 6- Gyro-Frequency Mop of the World Reference 35 



dependent varying inversely as the square of the operating frequency. At 130 
MHz ionospheric refraction is significant. The atmosphere influences Minitrack 
measurements during aircraft calibration and during spacecraft tracking. 


AIRCRAFT CALIBRATION 

1. Only the troposphere is a consideration for both radio and optical (star 
background and flashing light) refraction. 

2. The National Bureau >f Standards exponential atmosphere coupled with 
"flat earth" geometry is satisfactory for modeling aircraft calibration 
refraction 

3. For this model only the scalar quantity, surface refractivity, (optical and 
radio) is required. 


SPACECRAFT TRACKING 

1. Only the ionosphere is a consideration since the troposphere is corrected 
for inherently by the Minitrack calculation. 

2. The Ionosphere can be modeled using a Chapman type profile to obtain 
estimates of ionospheric induced bias. 

3. Recent top-side satellite soundings suggest an electron density decay 
with altitude somewhat less rapid than predicted by the Chapman model. 
An exponential decay in electron density above the height of maximum 
density appears to more closely approximate the true profile. 

4. The parameters - Maximum electron density, height of maximum density, 
am? scale height are required. 

5. The extent to which corrections are successful is highly dependent on 
validity of estimates for maximum electron density. 

6. Published prediction for maximum density used at the beginning and 
end of a given mor.Ci can be such + hat as much as 50% of ionosphere 
residual bias remains after corrections are applied. 
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